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Abstract. In the last century the non-perturbative regularization of chiral fermions was a long- 
standing problem. We review how this problem was finally overcome by the formulation of a 

i^ . modified but exact form of chiral symmetry on the lattice. This also provides a sound definition 

of the topological charge of lattice gauge configurations. We illustrate a variety of applications to 
^ \ QCD in the p-, the £- and the 5-regime, where simulation results can now be related to Random 

Matrix Theory and Chiral Perturbation Theory. The latter contains Low Energy Constants as free 

p^ ' parameters, and we comment on their evaluation from first principles of QCD. 
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CHIRAL SYMMETRY 



Chiral Perturbation Theory. Fermion fields can be decomposed into a left- and a 
>• \ right-handed component by means of the chiral projectors, ^l,r = 5(1='= Ys)^ ^ ^l,r = 

^yQ I ^j{^ TYs) ■ In massless (bilinear) theories these two spinor components decouple. In 

CN ■ particular the QCD Lagrangian at zero quark masses takes the structure 

O 

o 
o 



^ : ^QCD = ^L£>^L + ^i?^^i? + ^gauge , (1) 



where D is the Dirac operator, and the quark fields ^, ^ capture the Nf flavors and 
3 colors involved. This -Sfgco is invariant under global U{Nf) transformations of the 
quark fields, which can be performed independently in the left- and right-handed sector. 
The two complex phases represent baryon number conservation and an axial symmetry, 
cd I which breaks under quantization (axial anomaly). One assumes the remaining chiral 

flavor symmetry to break spontaneously, reducing the symmetry to a unbroken group of 
simultaneous transformations in both sectors, 

SUiNf)L ® SUiNf)R ^ SU{Nf)L+R . (2) 

Chiral Perturbation Theory (;tPT) deals with a field in the corresponding coset space, 
U G SU{Nf), which dominates the low energy behavior. As we add small quark masses 
iriq to -^QCD (which is allowed, since QCD is a vector theory), U represents NJ — l light 
quasi-Nambu-Goldstone bosons, which are identified with the lightest mesons involved. 
For simplicity we consider only two (degenerate) flavors, u and d, so that the field U 



represents the pion triplet.' ;f PT now uses an effective Lagrangian of the form [2] 

=^eff [t/] = ^Trid^U'^d^U] - ^Tr[[/ + 1/ V • • • (3) 

where the dots represent terms with more derivatives and/or higher powers of the explicit 
symmetry breaking parameter niq. All terms, which are compatible with the symmetries, 
are put into an energetic hierarchy; eq. (3) displays the leading terms. Each term comes 
with a coefficient, which is denoted as a Low Energy Constant (LEC), such as the pion 
decay constant F^^ (which can be measured experimentally) and the chiral condensate E 
(the order parameter of chiral symmetry breaking). 

As we have seen, ;^PT does have a direct link to QCD, hence it describes low energy 
hadron physics in a way manifestly related to the fundamental theory (in contrast to 
many other effective approaches). However, the LECs are free parameters in ;^PT; in 
this sense the low energy description is incomplete. If we manage to determine LECs 
directly from QCD, we obtain a more complete low energy theory. This is obviously 
a non-perturbative task, and therefore a challenge for lattice QCD: its simulation in 
Euclidean space is the only method to tackle QCD (and other quantum fields theories) 
systematically beyond perturbation theory. 

Lattice fermions. The lattice discretization of the gluon fields is conceptually un- 
problematic: the gauge action can be expressed in terms of small Wilson loops in a gauge 
invariant way. It has been a longstanding issue, however, to formulate lattice fermions 
such that they keep track of (approximate) chiral symmetry. For one flavor, the standard 
chirality condition is given by the anti-commutator {D, 75} = 0. The "naive" discretiza- 
tion of the Dirac operator yields for a free, massless fermion in momentum space the 
form D-a{p) = iy^sinp^ (in lattice units, i.e. if we set the lattice spacing a = 1). It is 
chirally symmetric, but it gives rise to 2^^ — 1 artificial poles of the propagator (inside the 
first Brillouin zone), in addition to the physical one at p = (ind dimensional Euclidean 
space). The Nielsen-Ninomiya Theorem states essentially that chirality and locality in- 
evitably entail fermion doublers, which would distort the result in lattice studies [3]. 
(Here "locality" means that the coupling between ^^ and ^y falls off at least exponen- 
tially in |.x; — jI; this assures a safe continuum limit). 

K. Wilson subtracted a discrete Laplacian A to construct the Wilson Dirac operator 
[4], Dw = D^ — jA. It is still local, and it sends the doubler masses to the cutoff scale, 
as desired. However, the additional term breaks chiral symmetry explicitly. Under gauge 
interaction it leads to (highly undesired) additive mass renormalization. Thus the chiral 
limit can only be approximated by a tedious fine-tuning of a negative bare quark mass. 

Conceptual progress was achieved at the end of the last century by deviating from the 
continuum form of chiral symmetry in a specifically harmless way: instead of inserting 
a local term for {D, 75} (as Wilson did), one now does so for {D~^, 75} — this does not 



' In this case, the symmetry breaking pattern (2) is locally isomorphic to 0(4) -^ 0{3). This property is 
very specific: among all conceivable types of chiral symmetry breaking only very few can be expressed 
by orthogonal Lie groups [1]. 



shift the poles in the propagator — e.g. by setting 

{D-\r5} = r5 ^ {D,r5} = Dr5D. (4) 

This is now known as the (simplest form of the) Ginsparg-Wilson Relation (GWR). Lat- 
tice Dirac operators with this structure are generated by block spinor Renormalization 
Group Transformations as the blocking factor diverges (perfect fermion action) [5, 6, 7], 
or by Domain Wall Fermions [8] (with the chiral modes attached to two domain walls, 
which are pulled apart in an extra "dimension"). By integrating out this extra "dimen- 
sion" one obtains the overlap fermion [9]. The latter can easily be re-derived from the 
GWR (4): assume D to be 75-Hermitian, D^ = YsDfs (which holds e.g. for Dw), so 
that the GWR can be written as D + D^ = D^D. This corresponds to the condition that 
A := D — 1 be unitary. For Aw = ^^w — 1 this is not the case, but we can enforce it by 

the transformation , 

Aov = Aw/ y A^Aw , £>ov = Aov + 1 , (5) 

which yields the "overlap Dirac operator" Dqv as a solution to the GWR [10]. Its 
spectrum is located on a circle in the complex plane, |A — 1| = 1, which reveals the 
absence of additive mass renormalization. (This also holds for the generalization to 
A = D — p, |A— p| =p,p>l, which has practical advantages in the interacting case.) 
The overlap operator is manifestly local as long as the gauge background is sufficiently 
smooth [11]. This property gets lost on very coarse lattices (a>0.17 fm), but if we 
simulate in the safe regime, a continuum extrapolation can be taken. 

Zero modes of a Ginsparg-Wilson Dirac operator are exact and they have a definite 
chirality. Hence the topological charge v of a gauge configuration can be defined [12] by 
adapting the Atiyah-Singer Index Theorem from the continuum. Thus v is defined as the 
difference between the number of zero modes with positive and negative chirality. (In 
actual Monte Carlo generated configurations only zero modes with one chirality occur.) 
In contrast to the formulation using Dw, operator mixing (on the regularized level) is 
under control [13]. This is very helpful for numerical measurements, for instance if one 
performs a fully non-perturbative Operator Product Expansion [14]. 

As a reason for these fantastic properties, M. Luscher pointed out that the Lagrangian 
is actually invariant under a lattice modified, chiral transformation with the infinitesimal 
form [15] 

1 1 

iI/^v|/(^l + eJl__Z)]yg) , vj/^(^i+£^5ji__£,jjij/ (toC>(e)). (6) 

In the continuum limit it turns into the standard chiral transformation, ^P — )■ ^(1 + £75), 
^ — )■ (1 + £75)^. However, the fermionic Lagrangian ^D^ is invariant to 0(e) even 
on the lattice, if the GWR holds. On the other hand, the functional measure &^! ^m 
is not invariant, which is exactly what it takes to reproduce the axial anomaly correctly 
[15, 16]. 

For applications also further properties — beyond chirality — matter. This motivates 
the substitution of the Wilson kernel Dw in the overlap formula (5) by an extended "hy- 
percube fermion", which is constructed from Renormalization Group Transformations 
and truncations [17]; this yields improved locality and scaling, as well as approximate 
rotation symmetry [18, 19]. In particular locality now persists on coarser lattices, which 
is profitable in studies of QCD at finite temperature [20]. 



APPLICATIONS TO TOPOLOGY AND P- , e- AND 5-REGIME 

We mentioned that the topological charge is well-defined when we deal with chiral 
lattice fennions (although all lattice gauge configurations can be continuously deformed 
into one another). This enables a sound numerical measurement of the topological sus- 
ceptibility Xt = i"^^) /y^ where v is the (aforementioned) topological charge, and V is 
the space-time volume. A high statistics study was presented in Refs. [21], which is very 
well compatible with our results [19], see Figure 1 (left). The continuum extrapolation 
amounts to Xt = (191(5) MeV)^. This supports the Witten-Veneziano scenario [22], 
which explains the heavy mass of the r] '-meson in part as a topological effect. This 
conjecture involves indeed the quenched value of ;^;. Latest direct measurements of m,^ 
and m,j/ were performed with 2-1-1 flavors of dynamical Domain Wall quarks [23]. 

XPT has been formulated in different regimes depending on the volume, which affects 
the counting rules for the energy hierarchy. 

p-regime. The standard setting, where finite size affects are small, is denoted as 
the p-regime [24]: L ^ l/niji (L is the 4d box length, and the inverse pion mass is 
the correlation length). Figure 1 (right) shows a measurement of F;,; that we performed 
(quenched) on a 12-^ x 24 lattice, where the gauge coupling was chosen such that the 
lattice spacing corresponds to a ^ 0.123 fm. The pion mass (and other hadron masses) 
can be measured from the exponential decay of correlation functions. Over a broad 
range, the measured value of F^i is clearly too high, but just at our lowest pion mass, 
m-,^ = 279(32) MeV, the right trend sets in, i.e. a decrease towards the value in Nature 
(F;i — 92.4 MeV). This calls for clarification with yet lighter pion masses, still closer 
to its phenomenological value of 135 MeV. However, at our lowest data point we 
already have Lnij,; ^ 2, hence at even smaller nij,; we are certainly outside the p-regime. 
Recovering it takes a much larger volume, and therefore much more computational 
effort. 

e-regime. As an alternative, we can simulate QCD in the e-regime [25], where 
Lmjt < 1- This setting is unphysical — experimentally we can't squeeze pions into 
such a tiny box. However, the finite size effects can be computed by ;tPT, and they are 
parameterized by LECs as they occur in infinite volume. Hence we can extract physical 
results nevertheless, avoiding the apparent quest for a huge volume [26]. In the £-regime 
the topological sectors play an essential role. 

Chiral Random Matrix Theory (RMT) provides a prediction for the density of the 
lowest non-zero Dirac eigenvalues A, (/ = 1 , 2, 3 . . .) in the £-regime [27]. More precisely 
it predicts the densities of the dimensionless variables zi = A/EV . If our data match the 
predicted shape, we can tune E for optimal agreement, and in this way evaluate E. Figure 
2 (left) shows the cumulative densities for zi as predicted by RMT in the sectors of 
topological charge | v | = 0, 1 and 2 (curves). Our data points, obtained in V ^ ( 1 .23 fm)"^, 
are in excellent agreement, if we insert E ~ (253 MeV)^ [28]. Also this result was 
obtained in the quenched approximation (which neglects sea quark contributions), but 
our result showed for the first time that this method to measure E is in fact applicable 
with chiral fermions in various topological sectors — see also Refs. [29, 19]. Recent 
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FIGURE 1. Left: The topological susceptibility of quenched QCD, measured with the index of chiral 
quarks. Our results [19] are consistent with the continuum limit of Refs. [21], which supports the Witten- 
Veneziano scenario for m„i . Right: F„ measured by two methods (straight pseudoscalar density correlator, 



and subtraction of scalar density correlator) at different pion masses in the p-regime (at a : 
For our lightest pion, niji — 279(32) MeV, a trend towards the value in Nature sets in [19]. 
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FIGURE 2. Left: The cumulative density of the dimensionless variable zi = AiEV, where Ai is the 
leading non-zero Dirac eigenvalue. We compare our data for Dqv eigenvalues in the e-regime (mapped 
stereographically onto IR+) to RMT predictions (lines) in the sectors with topological charge | v| = 0, 1 
and 2. We obtain very good agreement if we insert Z = (253 MeV)^ [28]. Right: The residual pion mass 
in the 5-regime as a function of the spatial box size L^ . The chiral extrapolation of our measured pion 
masses (with dynamical Wilson quarks) [40] follows closely the theoretical prediction of ;|^PT [38]. 

Studies with dynamical quarks [30] obtain with the same method a very similar value, 
E = (251(7) ]V[eV)3 (renormalized in the MS scheme at 2 GeV) . 

However, simulations of dynamical overlap fermions are not only computationally 
very expensive (the inverse square root in eq. (5) has to be computed by polynomials up 
to degree 0(100) in order to attain chirality close to machine precision), but they also 
face conceptual problems: the standard algorithm for dynamical quarks ("Hybrid IMonte 
Carlo") changes the topological sector only very rarely, so that direct measurements of 
full observables are difficult. IVIeasurements can be performed in fixed topological sec- 
tors. A method to derive from them an approximate result for the physical value (prop- 
erly summed over all sectors) has been suggested in Ref. [31] and tested successfully in 
the 2-flavor Schwinger model [32]. The tremendous efforts to simulate QCD with dy- 
namical overlap quarks are reviewed in Ref. [33]. In particular Xt is hard to measure in 
this case; for an indirect method we refer to Refs. [34]. 



We add that also Fji can be evaluated in the e-regime, in particular by matching 
measured correlators [35, 19] to ;^PT predictions, or by considering only their zero- 
mode contributions [36, 19]. 

5-regime. Let us finally mention yet a third regime where ;fPT has been worked 
out, namely the 5-regime [37]. Here the Euclidean time extent is long, but the 3d spatial 
volume, say L^, is small (L^ < l/m^). This prevents spontaneous symmetry breaking, 
hence even at vanishing quark mass — where we can refer to the current quark mass 
measured through the PCAC relation — the pion mass remains finite. The formula for 
the residual pion mass nf^^{Ls) in the chiral limit has been computed recently to next- 
to-next-to-leading (NNL) order [38].^ 

We performed p-regime measurements of m-n in spatial volumes in the range Ls ~ 
(1.6... 3.0) fm, and extrapolated the pion masses to the chiral limit [40]. The results 
are in remarkably good agreement with the formula for m^^(Ls). The latter involves F^i 
again, so in principle this is yet another way to measure a physical LEC in a unphysical 
regime (although in our case we already used F^i for the extrapolation). 

Moreover the NNL order also involves sub-leading LECs of ;fPT (coefficients to 
terms symbolized with dots in eq. (3)). Hence from precision results for nf^^{Ls) in 
the 5-regime one could determine even sub-leading LECs. That is useful in particular 
for the LEC denoted as l^ [40], the value of which is quite uncertain [41]. 

CONCLUSIONS 

Chiral fermions can be regularized on the lattice such that they obey a lattice modified 
version of chiral symmetry. Therefore they are now well-defined and tractable non- 
perturbatively, at least in vector theories. Thus the existence of light quarks — with 
masses far below Aqcd — is not that mysterious anymore (for reviews, see e.g. Refs. 
[42]). This formulation also provides a sound definition of a topological charge, which 
enables a neat measurement of the topological susceptibility in quenched QCD. The 
results support the Witten-Veneziano conjecture about the r^'-mass. 

We sketched applications of chiral lattice fermions in the p-regime, where finite size 
effects are small. Here we can measure for instance the light hadron spectrum as well as 
the PCAC quark mass, and Fti — one of the leading LECs in the ;t^PT Lagrangian. The 
LEC determination from the underlying theory (QCD) improves the status of ;^PT as a 
description of the low energy hadronic world in a way linked to first principles. 

In the e-regime we discussed the measurement of E — the other leading LEC — by 
relating the microscopic Dirac spectrum to Random Matrix Theory. Also here, and in 
the 5-regime, F^i can be measured, i.e. it is possible to obtain physical results even from 
unphysically small volumes. Moreover the measurement of the residual pion mass in the 
5-regime even has the potential to determine sub-leading LECs. 

Acknowledgments : I thank my collaborators in the works summarized here, and the 
organizers of the pleasant workshop in Mazatlan. 



This is very different from the p-regime, where finite size effects are suppressed exponentially [39]. 
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